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We compute the total cross section and the transverse momentum distribution for single charged 
slepton and sneutrino production at hadronic colliders including NLO supersymmetric and non- 
supersymmetric QCD corrections. The supersymmetric QCD corrections can be substantial. We 
also resum the gluon transverse momentum distribution and compare our results with two Monte 
Carlo generators. We compute branching ratios of the supersymmetric decays of the slepton and 
determine event rates for the like-sign dimuon final state at the Tevatron and at the LHC. 



I. INTRODUCTION 

Supersymmetry [l[ is the most promising and wideh 
studied solution to the hierarchy problem 0, H,_P 
of the Standard Model of particle physics (SM) 
It must be broken with a mass scale < 0(1 — 10 TcV) 
[1, M, El and thus should be testable at the Teva- 
tron and the LHC. When extending the SM by su- 
persymmetry, we must introduce an additional Higgs 
doublet, as well as double the spectrum. The most 
general renormalizable superpotential consistent with 
this minimal particle content as well as the gauge sym- 
metry of the SM is [H 03] 



W = W Pe + W Pe 
W P6 = hfiLiHjEi 



hfjQiH d Dj 
+ hfjQiH u Uj + [iRdHu , 

Wp B = XijkLiLjEk + !.,(), I), 
+ \'l jk U i D D k + n i L i H u . 



(LI) 



Here wc use the standard notation of [X 31 ] . in par- 
ticular: k = 1,2,3 are generation indices. This 
superpotential leads to rapid proton decay via the 
LQD and UDD operators [3| in disagreement with 
the lower experimental bound [Ta |. The supersym- 
metric SM thus requires an additional symmetry to 
stabilize the proton. It was recently shown, that there 
are only three such discrete symmetries which are con- 
sistent with an underlying anomaly- free U(l) gauge 
theory [lj| [l?], UK]: R-parity (or equivalently mat- 
ter parity), baryon-triality, and proton- hexality. It- 
parity prohibits Wp 6 and thus apparently stabilizes 
the proton. However, it allows dangerous dimension- 
five proton decay operators such as QQQL This 



* dreiner@th.physik.uni-bonn.de 
i sgrab@th.physik.uni-bonn.de 
-'mkraemcr@physik.rwth-aachcn.de 

SSince July 2006: MPI fiir Physik, Fohring er Ring 6, 80805 
Miinchen, Germany; trcnkcl@mppmu.mpg.de 



long-standing problem in supersymmetry is resolved 
by proton-hexality, Pq [T(|j a discrete Ze-symmetry, 
which acts like R-parity on Eq. (jl.ip but in addition 
disallows the dangerous dimension-five operators. 
An equally well motivated solution to the proton decay 
problem is baryon-triality, a discrete ^-symmetry, 
which prohibits the UDD operator in Eq. (|I.1[) . The 
baryon-triality collider phenomenology has three main 
distinguishing features compared to the Pg-MSSM 
[HI HU (minimal supersymmetric standard model) . 

1. The lightest supersymmetric particle (LSP) is no 
longer stable. 

2. Supersymmetric particles can be produced 
singly, on resonance. 

3. Lcpton flavour and lepton number are violated. 

These lead to dramatically different signatures at 
hadron colliders [2l|, [22j . It is the purpose of this pa- 
per to investigate specific baryon-triality processes at 
the Tevatron and the LHC in detail. We focus on 
the resonant production of charged slcptons and sneu- 
trinos via the operator LQD of Eq. (|L. L|) and their 
subsequent decay into observable final states at the 
Tevatron and the LHC. This is the only novel produc- 
tion process and it has the highest kinematic reach in 
the supersymmetric masses. 

Since the LSP is not stable, it need not be the light- 
est neutralino as in the MSSM, e.g. it could be a 
scalar tau p3| . However, the conventional supersym- 
metric parameter points considered in collider studies 
to date, the SPS points (24|, imply a neutralino LSP 
and we restrict ourselves to this case here. We there- 
fore consider the following reactions 



Ujd k -> q 
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Here ui, l~[ denote the sneutrino and the charged slep- 
ton, and xl, xt denote the neutralino and chargino 
mass eigenstates, respectively. The decays (|I.2b . d) are 
via the LQ-D-operator. The decays p.2b . c, e, /) are 
cascade decays typically via a gauge coupling. 
In 1997 ZEUS and HI observed a slight anomaly in 
their high-Q 2 data [IE Hl[, which could potentially 
be interpreted in terms of resonant squark production 
via the LQD operator [IE HE HI- However, after 
computing the NLO QCD corrections to squark pro- 
duction at HERA [3(| and also comparing with NLO 
QCD corrections to leptoquark pair production at the 
Tevatron [3lj ] . this interpretation was not consistent 
with the Tevatron data. We take this as a motivation 
to compute the NLO QCD and supersymmetric QCD 
corrections to the resonant slepton total and differen- 
tial production cross section at the Tevatron and the 
LHC. 

Resonant slepton production at hadron colliders via 
the LQD operator was first investigated in [32|, [33j . 
using tree-level production cross sections and the de- 
cays of Eq. (|T2l). f). In 0, HE HI] a detailed phe- 
nomenological analysis of E q. (II. 2b) was performed in- 
cluding a Herwig [13, HE l39l |4Q|| simulation, again 
with tree-level cross sections. The NLO QCD cor- 
rections to the cross section were first computed in 
[ill , EH . We go beyond this work in several aspects. 
First, we complete the calculation by including the 
SUSY-QCD corrections. These involve virtual gluinos 
as well as tri-linear scalar couplings due to the su- 
persymmetry breaking A-terms [43|. As we shall see, 
these can be substantial in certain regions of param- 
eter space and modify the QCD prediction by up to 
35%. We have also resummed the pt distribution of 
the final state gluons, c.f. EH. In order to facilitate 
the implementation of our results in experimental in- 
vestigations we present a detailed comparison with the 
Herwig and Susygen [|1] Monte Carlo programs. 
Furthermore, we present a detailed phcnomcnological 
analysis of the final state including the relevant decay 
branching ratios. Wc also present event rates at the 
LHC and Tevatron for the promising like-sign dilcpton 
signature using our improved calculation. 
In related work, resonant squark production at hadron 
colliders via the UDD operator at tree level was dis- 
cussed in [33| . The NLO corrections and the pr resum- 
mation has been studied in [45|. At HERA resonant 
squark production has been considered in [IE EE E3 • 
At e + e~ colliders one can have resonant sneutrino pro- 
duction via the L i Lj E-\ -operator. This has been con- 
sidered in [II EE El. 

At hadron colliders resonant production proceeds via 
the operators \' i ^ k L i QjDi : . However, for j, k = 2,3 
the production rate is suppressed due to the lower lu- 
minosity for sea quarks. The current best bounds for 
the first generation quarks are given by [IE HE HU . 



Coupling 


\' \' \' 

A lll A 211 A 311 


best bound 


0.0005 0.059 0.11 



The last two bounds scale with the relevant scalar 
fermion mass (mj/lOOGeV). The first strict bound 
arises from ncutrinolcss double beta decay searches 
and has a more complicated rrij dependence. 
The outline of our paper is as follows. In Sect. [Hi we 
compute the NLO resonant slepton production cross 
section at hadron colliders, including both virtual and 
real diagrams. We also employ this calculation to de- 
termine the running of the parameter Ay fe . We then 
first compute the partonic cross section and then nu- 
merically the hadronic cross section. In Sect. IIII1 we 
investigate the transverse momentum distribution of 
the produced sleptons at NLO. In particular we per- 
form a pt resummation of the gluon distribution and 
then compare with the Herwig and Susygen Monte 
Carlo generators. In Sect. II VI we study the branching 
ratios of the various slepton decays and determine the 
resulting like-sign dimuon event rates at the LHC and 
the Tevatron. In Sect. [V] we summarize our results 
and conclude. 



II. TOTAL PRODUCTION CROSS SECTION 

A. LO Processes 

We first consider the single slepton production pro- 
cesses at leading order, Fig. [TJa), 



dk{p) + dj/uj(p') -> i>i/£u(q) 
dkip^ + dj/u^^ut/iUq). 



(II.1) 



The relevant parts of the Lagrangian for the produc- 
tion at hadronic colliders follow from the superpo- 
tential Eq. (|I.1[) . Expressed in terms of the (four-) 
component fields and using the projection operators 



P, 



R/L 



i(l ± 7 5 ) we have 



C-LQD ^>K 3 k {&i d k Ph dj - hi d k P L Uj) 

+ \% k (DtdjPRdk-I+^Pndk). 



(II.2) 



We obtain for the total partonic single slepton 
(charged slepton or sneutrino) production cross sec- 
tion (32|, [33[ 



'LO 



\K jk \^ 

12s 



8(1 - r) = o-o 8(1 - t) (II.3) 



with r = m 2 /s, where rh is the slepton mass and 
the partonic center-of-mass energy. 



B. NLO SUSY-QCD Corrections 

The cross section at Q(a s ) comprises radiative 
corrections to the quark-antiquark process includ- 
ing virtual gluon, gluino, quark and squark ex- 
change, Fig. [TJb), (c), (d), (e), real gluon radiation, 
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Fig. QJf), and Compton-like gluon initiated subpro- 
cesses, Fig. [ljg). In the following, we shall concen- 
trate on positively charged slepton production. The 
production of sneutrinos can be treated analogously. 
Throughout this section, generation indices are sup- 
pressed. 

We use dimensional regularization to isolate ultravio- 
let, infrared and collinear divergences, and adopt the 



MS scheme for the renormalization of the coupling 
A' and the factorization of the collinear divergences 
through a redefinition of the parton distribution func- 
tions. 

The calculation of the NLO QCD corrections from 
virtual gluon exchange and real gluon emission is 
straightforward. For the partonic qq' cross section at 
NLO we find 



-NLO. QCD _ 
a qq' 



a s „ /3 , U % 3 , (il 7T 2 



cr — Of 

7T 



1 



ln^ + 2(l + r 2 ) 
T A^f 



ln(l - r) 



(1-t) 



(II.4) 



and for the quark-gluon scattering contribution 



-NLO -NLO 

a ig =a <i' 9 



+ [r 2 + (1 - r f 



' m 2 , (1 — t) 
In -5- + In 

. Mf t 



(II.5) 

where Cf = 4/3 and r = m 2 /s. The couplings A'(/zf) 
and ol s (hr) are evaluated at the renormalization scale 
fin, and /^f denotes the factorization scale. Our re- 
sults agree with those presented in [4l|, E3] . 
For a complete cross section prediction at 0(a s ) one 
has to include SUSY-QCD corrections through virtual 
gluino and squark exchange. These corrections have 
not been considered in previous calculations. The 
SUSY-QCD corrections at NLO arise from the ver- 
tex correction, Fig. O^d), and the quark self energy, 
Fig. HJe), and are infrared- finite. 

The interaction of quarks q, (right-) left-handed 
squarks qi, (c[r) and gluinos g is described by [52l. |53| 

C qqg = -V2g s t a bd J2 P l Qf 9Lf + Qf PR ~9a f L / 

f=u,d 



g a Pr q s q Rf - q s Pl g a q Rf 



(II.6) 



Here, / denotes the quark flavor, g s the QCD coupling, 
b, d (a) arc the indices of the fundamental (adjoint) 
representation of SU(3)c, and t^ d are the generators 
of the group. Squark mixing between left- and right- 
handed squarks is neglected since only light quarks 
in the initial state contribute to the hadron collider 
cross section. Moreover, mixing effects between squark 
generations are excluded to avoid FCNC processes. 



slepton are needed. They can be extracted from the 
P fi -violating parts of the general soft SUSY-breaking 
Lagrangian [20j 



£ ST D ~ X 'tjk A (pi dhj d* Rk ~ hi u L j d*R k 



\ /* A * 

A ijk A 



d* L j dRk 



~Li U*Lj dRk 



(II.7) 



We have assumed here that the soft breaking terms 
have a universal dimensionful parameter A and are 
proportional to the dimensionless coupling constant 
Ay ft of the superpotential [Hj]. 

Due to angular momentum conservation only those di- 
agrams contribute where the incoming quark and anti- 
quark have the same helicity. Therefore, for a specific 
slepton production process, only a single vertex correc- 
tion diagram contributes at NLO which is described by 
Eq. d^- 
Using standard notation [5|| for scalar loop integrals, 
the SUSY vertex contribution to the NLO cross section 
can be written as 



-NLO 
"SUSY, 



•LO "a 



C F Am s Co(-p,p',m g ,mg L ,mq R ) , 



(Hi 



where p,p' are the momenta of the incoming quark 
and antiquark, m g is the gluino mass, and rriq L , m g - H 
are the masses of the left- and right-handed squark 
and antisquark running in the loop. 



1. Vertex correction 

For the SUSY-QCD vertex corrections, also the in- 
teraction vertices of two squarks and a (left-handed) 



2. Quark Self Energy 

SUSY-QCD self energy contributions arise from di- 
agrams of the form Fig. QJe). Since the SUSY part- 
ners of left- and right-handed quarks can have different 
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FIG. 1: Generic Feynman diagrams for resonant slepton 
production, (a) LO; (b) QCD vertex corrections; (c) QCD 
self energy; (d) SUSY-QCD vertex correction; (e) SUSY- 
QCD self energy; (f) QCD real gluon radiation; (g) QCD 
Compton process. 



masses, the SUSY self energy corrections to the left- 
and the right-handed quark field differ in general. 
The field strength renormalization constant Z is 



Z- 1 = 1- |Ev-£vM . ( IL9 ) 



where the vector part of the self energy £y can be 
written in terms of the coefficient function B\ [55l as 



A' 



2tt 



C F Bi(-p,m.s,m,- 1/B ). 



(11.10) 



The UV pole of the integral is, in the MS scheme, 

1 



Bi\uv = -(4^) e r(l + e) 



2e' 



(11.11) 



where the number of space-time dimensions is TV = 
4-2e. 

The SUSY-QCD self energy contributions to the NLO 



production cross section follow directly, 

» NLO _ -LO a s r < 
Z7T 

x (^Bi(~p,mg,m,q L ) + Bi(-p,mg,mg R ) - 2Bi\uv\ 

(11.12) 

This completes the computation of resonant slepton 
production at NLO. 

C. The running coupling A' 

In this section, we briefly discuss the 0(a s ) running of 
the LQD coupling A' including both QCD and SUSY- 
QCD corrections. 

The bare LQD coupling A' is related to the renormal- 
ized coupling via 



A 



zav 



(11.13) 



where Z is the renormalization constant for the cou- 
pling A' and the scale \x is introduced to keep the renor- 
malized coupling dimcnsionless. Note that the only 
scale dependence of Z is carried by a s and, therefore, 

dZ _ dZ da s _ f3(a s ) dZ (II 14) 

d/i da s d\i [i da s ' 

where the /3-function is defined as /3(a s ) = —\xda s j d\x 
with f3(a s ) = 2ea s + 0{a 2 s ). We thus obtain for the 
/3-function of the LQD coupling A' 



/3(A') 



dX' 
dfx 



A'/?K)^ + A'e. 
Z da s 



(11.15) 



Considering only QCD corrections, the renormali- 
zation constant Z depends on the (identical) renorma- 
lization constants for the lcft-/right-handcd quark field 



7j l IR an d on the vertex renormalization Z\* . They are 
Z x , = 1- ^ C F (4^) e r(l + e)i, (11.16) 



J L/R 



^C F (4^) e r(l + e)i, 
7T e 

a 1 

i-^tM^ra + e)-. 



(11.17) 



Taking also the SUSY-QCD corrections into account, 
SUSY self energy contributions Z L j R alter the renor- 
malization constant 

Zl/r = l-^C F (47r) £ r(l + e)i. (11.18) 
47T e 

Explicitely, including only QCD corrections, 

Z = ^L 
Z L 

and including QCD and SUSY-QCD corrections 



1 



^CV(47rrr(l + e)i, (11.19) 
47T e 



Z 



Z L Z L 



1- ^CV(47r) e r(l + e)l. (11.20) 

Z7T e 
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For the /3-function Eq. 15[) one obtains 

/3(A') = (/3 A Q , CD + /?f, USY ) A'g + 0(X'all (11.21) 

with the coefficients /^, CD , /3fF SY related to QCD and 
SUSY-QCD effects, respectively, 



„QCD or< oSUSY 
QCD 



(11.22) 



such that in total f3 tot 



By integrating Eq. (|II.15|) and inserting the definition 
Eq. (|II.21|) , the NLO running of the LQD coupling is 
obtained as 



2C F 



A' (/xo) 



Mo 



(11.23) 



in agreement with the results presented in [23j, |56| ; see 
however [931. 



D. Decoupling of heavy SUSY particles 




m [TeV] 



FIG. 2: Single slepton production cross section at the Teva- 
tron (lower curves) and at the LHC (upper curves). The 
LQL>-coupling is set to A$ u = 0.01. The CTEQ6 LO/NLO 
PDFs have been used, and renormalization and factoriza- 
tion scales have been identified with the slepton mass fa. 



In the limit of heavy SUSY particles rriq L , m, H , m§ = 
M and M — ► oo, the NLO SUSY-QCD vertex correc- 
tions vanish 

&suin, vc = -^lo — Cp AM C Q (-p, p', M, M, M) 



0, 



(11.24) 



while the self energy contributions are logarithmically 
enhanced 

*susy, so = *lo g C F [B! (-p, M, M) + B 1 (-p, M, M)] 



_cr LO — t-F In 

u 2tt M 2 



(11.25) 



The large logarithms can be absorbed by a redefinition 
of the coupling A' according to 



A' = A' (mo) 
1 



l-^/^ln^ 
Air Mo 



A' + AA' = A'(mo) [l - ^ fr? In 4 - ? ^ ^ ^ 

L 47T ^5 47T M 2 

(11.26) 

Then, the running of A' is determined only by the QCD 
corrections (light particles) 



d{\' + AAQ 
d In /i 2 



-A'(W))g[^ + CF] 
-A'(Mo)^ CD 



(11.27) 



Expanding A' = (A' + AA') - AA' in <t lo in Eq. (|lL25)) 
cancels the term oc ln(// 2 /M 2 ) so that in the limit 
M — * oo the SUSY-QCD corrections vanish. Note 
that we will not adopt the decoupling scheme de- 
scribed above in the subsequent numerical analysis as 
we do not consider scenarios where squark and gluino 
masses are much larger than the slepton mass. 



E. Total Hadronic Cross Sections 

In this section, we will present numerical results for 
the total hadronic cross section for resonant charged 
slepton production at the Tevatron and at the LHC. 
We assume only a single non-zero coupling = X' iud 
and set this coupling A' = 0.01 at all scales. The 
cross sections are directly proportional to |A'| 2 such 
that results for other values of the coupling are easily 
obtained by rescaling. For the calculation of the PP 
and PP cross sections we have adopted the CTEQ6L1 
and CTEQ6M (5?| parton distribution functions at 
LO and NLO, corresponding to A LO = 165 MeV and 
Af s = 226 MeV at the one- and two-loop level of the 
strong coupling a s (fj,u), respectively. The renormali- 
zation and factorization scales have been identified 
with the slepton mass, fiR = fip = m. 
The hadronic production cross section for the process 
ud — > if at the Tevatron and at the LHC is shown as 
a function of the slepton mass in Fig. [2] We compare 
the LO result with the NLO results with only QCD 
and both QCD and SUSY-QCD corrections included. 
The NLO cross section for producing sleptons with 
a mass m = 100 GeV is approximately 1 pb at the 
Tevatron and 10 pb at the LHC, decreasing rapidly 
with increasing slepton mass. The slepton discovery 
reach at the Tevatron is thus limited by the compara- 
bly small cross section, with only one expected event 
at m = 600 GeV (for A- n = 0.01) for an integrated 
luminosity of JdtC = lfb -1 . The cross section is 
more favourable at the LHC, where a 600 GeV slep- 
ton would be produced with a ss 50 fb, leading to 
500 events per year even for the low luminosity phase 
with C = 10 fb" 1 /year. 

To demonstrate the impact of the SUSY-QCD cor- 
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FIG. 3: SUSY-QCD A-factor for resonant slepton production via A-u = 0.01 at the LHC, y/S = 14 TeV, for three sets 
of squark and gluino masses (rriq L = rriq R = m g ). Different values for the SUSY breaking parameter A are chosen as 
indicated. For comparison we give also the QCD A-factor (solid lines). The CTEQ6 LO/NLO PDFs have been used, 
and renormalization and factorization scales have been identified with the slepton mass m. 
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FIG. 4: The same as Fig. |3]but for the Tevatron. 



rections more clearly, we consider the i^-factor K = 
cnlo/clo denned as the ratio of NLO to LO cross 
section, with all quantities calculated consistently in 
lowest and next-to-leading order. In Fig. [3] ([4]) we dis- 
play the QCD as well as the SUSY-QCD X-factor at 
the LHC (Tevatron) for three different sets of squark 
and gluino masses and for three different values of the 
SUSY breaking parameter A. 

The QCD A'-factor shows an enhancement of the LO 
cross section of up to 35% at the LHC and of up to 
50% at the Tevatron. For scenarios where the trilincar 
SUSY breaking parameter A vanishes, the SUSY-QCD 
corrections are moderate. They enhance the cross sec- 
tion prediction by up to 10% for light sleptons and 



reduce it by 5 — 10% for heavier sleptons. For non- 
vanishing A there are additional contributions due to 
the SUSY-QCD vertex corrections which can lead to 
a significant further enhancement or reduction of the 
cross section by up to 15% for rriq.g = 300 GeV. For 
light squark masses, it can be as much as 35%. One 
observes the usual resonance peak when the slepton 
mass is twice the mass of the squarks in the vertex 
loop. 

We remark that the inclusion of the NLO-QCD cor- 
rections reduces the factorization scale dependence of 
the theoretical prediction as expected. At the Teva- 
tron, we consider the slepton mass interval from 100 
GeV to 1 TeV. A variation of the factorization scale 
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between fJ-F = 2m and /ip = m/2 then changes the 
LO cross section by up to ±20%. This uncertainty 
is reduced to less than ±10% at NLO. At the LHC, 
we consider slcpton masses between 100 GeV and 2 
TeV and again vary the factorization scale between 
fip = 2m and [if = m/2. The LO uncertainty is 
then up to ±10% which is reduced to less than ±5% 
in NLO. 

We conclude this section by tabulating some selected 
cross section predictions for slepton production with 
fa = 300 GeV at the Tevatron or the LHC via the ud 
initial state in Tablcfl] We show results for the de- 
fault CTEQ6M_parton density functions (PDFs) and 
for the MRST [SI, parametrization. We find that 
the difference between the two PDF parametrizations 
is less than approximately 5% for the NLO cross sec- 
tion predictions. 



simulations. 



A. Soft and Collinear Gluon Summation 

The summed differential cross section, valid in the re- 
gion of soft and/or collinear emitted partons, is given 
by M 



~ 2 roc 



dp T dy 



dbbJ (bp T )W{b) , (III.l) 



where S is the hadronic cm. energy, Jo the zeroth 
order Bessel function, and <jq = <To(s = fn 2 ) is the 
leading order cross section defined in Eq. (|II.3[) . The 
Sudakov-like form factor W(b) is given by 



III. SLEPTON TRANSVERSE MOMENTUM 
DISTRIBUTION 

In this section we compute the transverse momentum 
(j>t) distribution of the produced sleptons. At Born 
level, momentum conservation enforces the produced 
sleptons to have zero pr- At 0(a s ) the sleptons may 
have non-zero px due to an additional parton in the fi- 
nal state. The differential hadronic cross section in px 
at 0(a s ) is given in Appendix |A"1 However, for pt <C 
m, with fn the slepton mass, the emitted parton be- 
comes soft and/or collinear to one of the incoming par- 
tons leading to large logarithms in the differential cross 
section. The perturbation series in a s breaks down 
and has to be replaced by a series in a™ log™ (fn 2 /p T ) 
with m, n £ N and m = 0, . . . , 2n — 1. We shall employ 
the Collins, Soper, and Sterman (CSS) formalism to 
resum the large logarithms [6(| [HI, [62] ■ 
The pt spectrum of the sleptons has also been stud- 
ied in [42j. We present the main formulae and extend 
the studies of [42| by a comparison with Monte Carlo 



W(b) =cxp 

'(6o/b) 2 1 



dq 2 a s {q 2 ) 



2^ r 



2- 



K x In. 



f qA {xl)h>B{x%) + {xl~xl) 



, (III.2) 



where Xy 2 



& v \J fh 2 / S are the parton momentum 
fractions in the limit px/m — > and y denotes the 
slepton rapidity. The f q A, fq'B are the parton density 
functions for hadrons A, B, respectively. The coeffi- 
cients Ki, K2 are in general functions of a s . However, 
as we discuss below, we are here only interested in the 
exponential at 0(a s ), i.e. K12 are constant. The in- 
tegration in Eq. pil.2|) has to be performed over the 
impact parameter b, the Fourier conjugate of pt- We 
perform the integration analytically by using the LO 
running of a s . The PDF's are evaluated at a fixed fac- 
torization scale (J,p, where the most convenient choice 
[62l ] is (If = bo/b with 60 = 2e _7E and je the Euler- 
Mascheroni constant. 

We calculate the coefficients K12 by expandin g (|III.1I) 
up to 0(a s ). After some algebra one obtains [63l| 



^asym &q ^2 ^ 
7T 



K t ln-5- 

Pt 



f qA {xl)f q , B (xl) 



dpT dy S pt 

+ flA(4)(P?<°fiB)(4)+ J2 h'A{x 1 ){P qi of iB )(x° 2 ) 



i=q,g 



*=q,g 



(III.3) 



where Pg'i, Pqi are the Altarelli-Parisi splitting func- 
tions (64[ , and the convolution is defined as 

(/•»)W-jff /(•)«(?)■ tnu) 

We can then compare the soft- and collinear diver- 
gences of this asymptotic cross section with those of 



the perturbative result in the limit pt — > 0. We obtain 
A'i = 2C F and K 2 = -3C F with C F = 4/3. 
It was pointed out in [(55[ that W(b) is ill-defined when 
b > 1/Aqcd, because confinement sets in. We factor 
out the non-perturbative part by replacing W(b) by 



8 





CTEQ 6 


MRST 04 








SUSY-QCD 






SUSY-QCD 




LO 


QCD 


A = TeV 


A = 1 TeV 


LO 


QCD 


A = TeV 


A = 1 TeV 


a [fb] 


25.72 


38.43 


39.23 


40.22 


25.55 


37.73 


38.52 


39.50 


Te vat ron 


















a [fb] 


358.7 


467.0 


478.0 


491.6 


347.4 


471.2 


482.4 


496.2 


LHC 



















TABLE I: Numerical results for the production cross section for a slepton of mass m — 300 GeV at the Tevatron (y/S = 
1.96 TeV) and the LHC (y/S = 14 TeV). Slepton production via the ud initial state is considered, with A^n = 0.01. All 
scales are fixed to the central scale /i = in. S quark and gluino masses are chosen as m,q L = mg„ = m,g = 600 GeV. Results 
are presented for the CTEQ6L1 /CTEQ6M [57| and the MRST2002LO/MRST2004NLO [UlSl PDF parametrizations. 



G2. hi 



Results 



W(b) W{K) e - s ""^ , 

h h 
o* = — ; 

yfl + V/b. 



max 



(III.5) 
(III.6) 



which smoothly cuts-off the region b > 6 max , and our 
choice of 6 max is given below. S np (b) parameterizes 
the non-perturbative part and has to be determined 
by experiment. We follow the approach of (u7j . where 
g— S np (b) j g approximated by a Gaussian with adequate 
parameters fitted to Drell-Yan data resulting in 



S„ p (b) 



With 6 max 
finds gi 



9i + 92 In 



2Qc 



g 1 g 3 lii{100x 1 x 2 ) . 

(III.7) 



0.5 GeV" 1 and Q = 1.6 GeV, [67 



+0.05 
-0.04- 



0.21+°;°! GeV 2 , g 2 



0.68 



+ 0.01 
-0.02 



GeV z , g 3 



-0.60 

to factor out nonperturbative physics. 



Yang et al. [42| used the approach of [6£ 



B. Matching 

At low pt, the summed cross section describes the 
spectrum accurately compared to the perturbative 
cross section at C(a s ), Eq. (|A.2j) . All terms singular 
in px are included. For large pt, also the non-singular 
terms must be considered. They correspond to the dif- 
ference of the perturbative and the asymptotic cross 
section at (D(a s ). For large pr the summed (jlll.lj) and 
asymptotic (|III.3|) cross sections are negative, which is 
unphysical. Therefore, we suppress their contributions 
with an empirical matching function (69l . [70| : 

dV csum d 2 cr as >" 



d PT dy 1 + ( PT /p^ atch ) 4 



dp T dy 



dprdy 
dp T dy 

We choose p™ atch = rfi/6 at the Tevatron andp^ atch 
m/3 at the LHC. 



The full transverse momentum distribution for the 
case of a (left-handed) slepton at the Tevatron (LHC) 
is shown in Fig. O (Fig. [5]). The y-dcpcndcnce has 
been integrated out and the renormalization and fac- 
torization scales of the asymptotic and perturbative 
cross section have been set equal to fi — yjp\ + m 2 /2. 
Again we set X' in = 0.01 and use the CTEQ6M PDFs 
[57l |. We normalize the integral of the transverse mo- 
mentum distribution to the total hadronic cross sec- 
tion including QCD corrections. In the low-py region, 
the perturbative and the asymptotic cross section can- 
cel each other and the spectrum is described well by 
the summed prediction. At high py, i.e. pt > 50 GeV 
(100 GeV) at the Tevatron (LHC), the summed and 
asymptotic contribution start to fall off and the distri- 
bution is to a good approximation given by the pertur- 
bative cross section. The summed cross section peaks 
at p T = 3.1 GeV (p T = 3.4 GeV) at the Tevatron 
and at p T = 4.3 GeV (p T = 5.0 GeV) at the LHC 
for m = 200 GeV (m = 500 GeV) 
Yang et al. [42j derived the transverse momentum 
distribution for sneutrinos with masses of 200 GeV, 
400 GeV and 6 00 Ge V. They did not use a match- 
ing procedure (|III.8[) and did not normalize the in- 
tegral of the transverse momentum distribution to the 
total hadronic cross section including QCD correc- 
tions. To facilitate the comparison with [42| we have 
also calculated the pt distribution of the sneutrinos 
without a matching function and without normalizing 
to the total NLO cross section. In the phenomeno- 
logically relevant region, i.e. where pr<50GeV, 
there is good agreement with [Q. At high pt we 
find discrepancies. Our differential cross section at 
Pt = 100 GeV {p T = 200 GeV) is 75% (30%) larger 
(smaller) compared to [42| for a sneutrino with a mass 
of 200 GeV at the Tevatron (LHC). These discrepan- 
cies are not particularly relevant phenomenologically, 
because the differential cross section is between two 
and three orders of magnitude smaller at pt = 100 
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25 50 75 100 

p T [GeV] 



FIG. 5: Differential transverse momentum distribution of 
the slepton at the Tevatron for the ud initiated process via 
X'ni- The two sets of curves correspond to slepton masses 
of 200 and 500 GeV. The different lines of each set de- 
note the perturbative (dot dashed), asymptotic (dashed), 
summed (dotted) and full (solid) differential cross section. 



GeV (p T = 200 GeV) at the Tevatron (LHC) com- 
pared to the peak region. 

Within this computation, there are some sources of 
theoretical error. First, the choice of the scales \If 
and fiR and the empirical nature of the matching pro- 
cedure is ambiguous. These uncertainties mainly effect 
the high-px region, where there are only few events. 
Second, the parameters describing non-perturbativc 
physics, cf. Eq. (|III.7j) . are also not unique. Third, 
there are uncertainties from the choice of the PDFs. 
In order to quantify the latter, we have adopted the 
MRST2004NLO [53] PDF set for comparison. The 
maximal difference between the two pr distributions 
is at the peak. It is less than 0.9% (4.6%) at the 
Tevatron (LHC) for m = 200 GeV. Finally, there is 
of course an error due to higher order corrections of 
O{oi). 



D. Comparison with Monte Carlo Generators 

Next, we compare our analytic results for the slepton 
transverse momentum distribution with the distribu- 
tions predicted by two Monte Carlo event generators: 
Herwig 6.510 [37|, [H, [H Eland Susygen 3 Q, 
which uses the Pythia 6.205 |7ll . [72j parton shower. 




50 100 150 

PT [GeV] 

FIG. 6: Same as Fig.[5]but for the LHC. 



The Monte Carlo programs generate the pr distribu- 
tion from the LO process through the parton shower, 
for which we use the default settings. In Fig. [7] (Fig. [5]) 
we compare the transverse momentum distribution as 
predicted by the analytic resummation with Herwig 
and Susygen. All curves are normalized to the NLO 
QCD cross section. For both the Tevatron and the 
LHC, the spectrum generated with Susygen (through 
the Pythia parton shower) is significantly softer than 
that of Herwig. The analytic summation lies in be- 
tween the two Monte Carlo programs. The discrep- 
ancy between the Herwig and Pythia parton show- 
ers is a known problem, which also occurs, for exam- 
ple, for resonant single Z-boson production (73[. Our 
analytic results might help to improve future versions 
of Monte Carlo generators for resonant scalar particle 
production. 



IV. SLEPTON DECAY 

In order to determine the final state particles observed 
in a collider experiment as well as their event rates, 
we must discuss the possible decays of the resonantly 
produced (left-handed) sleptons. Since all SUSY par- 
ticles, even the LSP, are unstable in Pg-violating SUSY 
models, their decay signatures are crucial for experi- 
mental analyses. Here we shall assume a neutralino 
LSP. Note that in Pg-violating SUSY models one can 
also have a scalar tau LSP which dramatically alters 
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TABLE II: Relevant masses of the mSUGRA spectrum for the six Snowmass points, i.e. those which are kinematically 
accessible in the slepton decays. All masses are in GeV, rounded to the nearest whole number. 
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FIG. 7: Transverse momentum distribution of the slepton 
at the Tevatron generated with Herwig (solid histogram) 
and Susygen (dashed histogram) compared with the CSS 
formalism (solid line). The slepton mass is set to 200 GeV, 
and the curves are normalized to the NLO QCD cross sec- 
tion. 
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FIG. 8: Same as Fig. but for the LHC. 



the potential signatures [23|. 

Apart from the inverted production process, the slcp- 
tons can decay through gauge interactions. Neglect- 
ing mixing between left- and right-handed sleptons the 
dominant decays are, cf. Eq. (|I.2j) : 



djd k — > v* 



djdk, 

*-i Xn J 



Ujdk 



d k - ° + 



Ujd k , 

^1 Am ' 



(IV.l) 



In order to quantitatively discuss the possible signal 
rates, we consider the parameters of the six mSUGRA 
SPS points [24 L We calculate the spectra using the 
Spheno 2.2.3[Z3 online calculator [7l|, [zl, Jz| with 
the SM inputjz|: m t = 174.2 GeV, m b {m b ) m = 4.2 

GeV and af s (M z ) = 0.1176. The relevant masses 
of the SUSY particles are shown in Table [IT] in GeV. 
We have only listed those masses which are kinemat- 
ically accessible in the direct slepton decays of the 
given spectrum. In these scenarios the lightest neu- 
Possible decay modes via the 
k interaction are [H, [zl Hfl, 13 : 



tralino, is the LSP 
Kjk L iQj D 




Uidjdk 
Vidjdk 



(IV.2) 



The heavier neutralinos X2 3 4 an d the charginos X12 
may also decay via A^- fe but these channels are strongly 
suppressed compared to the gauge coupling decay 
modes into the LSP. For example, the branching ratio 

(BR) for the decay X2 — * [i + ud with X' 211 = 0.01 is 
less than 2 x 10 -6 in the five SPS scenarios, where the 
slepton can decay into an on-shell 
From the multiple decay modes we concentrate on 
those which lead to final states with two like-sign lep- 
tons; a signature with low background [36l . l82l . [83l . [8J, 
|85| . Here we assume A' 211 7^ 0, since dimuon final 



v 211 

states are of great interest for experimental analyses 
[S6l . H3, [H, For our numerical studies, we ne- 

glect decay modes with branching ratios BR < 0.5% 
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cessible. 

~+ _^ ^+ In the case of e and v e production, i.e. A' 1X1 7^ 0, one 

obtains the same number of events for a given value 
u + ud of A', however the bounds on A' m are much stronger 

[20i| . If \' 311 7^ there are more possible decay modes, 
fi + (i + X23 4, because left-right mixing in the stau sector can not 

be neglected. Both mass eigenstates n and f 2 can be 
produced via the L3Q1-D1 operator and the heavier 
u d mass eigenstate can decay into the lighter one via f 2 — * 

fiZ° and r 2 — > fih°, if kinematically allowed. This 
results in more decay chains compared to those for jl 
+ _^ ,-, -c+ and v Ll . In addition, the charge of the r-jcts has to be 

, _ reconstructed. 
*i At the LHC, the beam energy will be high enough to 

produce single sleptons at a high rate and an excess 



A' 



(up/j, )(u ud) 



.+ . 



of like-sign lepton pairs would be a hint of baryon- 



" ^- 1 ' 2 triality supersymmetry. Since the event rates scale 

W~ xi with |A- jfc | 2 for A- jfc much smaller than the gauge cou- 

A ' plings, we can roughly estimate that the magnitude of 

^ ^ ud A 211 can be tested down to 10~ 2 — 10~ 4 , depending on 

the scenario. 



TABLE III: Relevant decays of resonantly produced single 
smuons jl and smuon sneutrinos i> M into like-sign muon 
pairs of positive charge. 

and where there are less than 0.1 events for the given 
luminosity. The relevant decay chains leading to a a 
like-sign muon pair are shown in Table IIIII 
There are two distinct baryon triality scenarios. If 
X'211 is comparable in size to the gauge coupling, the jl 
and will have a significant decay width through the 
L2Q1D1 operator. In contrast, if A 211 is much smaller 
than the gauge couplings, the jl and decay through 

Pg-conserving channels. For example, the BR(jl + — > 
ud) is less than 1.1% for A 211 = 0.01 in all six scenarios. 
Baryon triality then manifests itself mainly through 
the LSP decays. Note that the BRs of the LSP at tree 
level are independent of A 211 . For the event rates, we 
choose A 211 = 0.05 and A' 211 = 0.01 in the following. 
In Table HV| we give the expected number of like-sign 
dimuon events (/z ± /x ± ) for the LHC (Tevatron) with an 
integrated luminosity of 10fb _1 (lib -1 ). We use the 
NLO cross sections including QCD and SUSY-QCD 
corrections with the tri-linear SUSY breaking param- 
eter A set to zero. Apart from the jl, and LSP 
decays, we use the BRs obtained by Spheno 2.2.3 

If decays into the heavier neutralinos X3 4 an d chargino 
X2 are kinematically forbidden (SPSla, lb, 3, 5), the 
jl decays dominantly into two like-sign muons associ- 
ated with two jets; contributions from the decay 
are negligible. Here all final state particles can be 
detected and the slepton and gaugmo masses can be 
reconstructed. If the and jl are heavier than all 
neutralinos and charginos (SPS2, 4), then the de- 
cay gives the main contribution and also decays with 
an on-shell Z° in the final state are kinematically ac- 



V. CONCLUSION 

We have argued that in preparation for the LHC, the 
baryon-triality interactions must be considered in de- 
tail. We have computed the full NLO corrections to 
resonant slepton production at hadron colliders. We 
have found significant changes in the total and dif- 
ferential cross sections. In particular, we have com- 
puted the SUSY-QCD corrections for the first time 
and found that the effects can be large in specific re- 
gions of parameter space, modifying the QCD predic- 
tion by up to 35%. The factorization scale dependence 
of the cross section at the LHC (Tevatron) is also sig- 
nificantly reduced to about ±5% (±10%) at NLO. We 
have used our computation to determine the running 
of the A'-coupling without and with the decoupling of 
heavy SUSY particles. 

Next, we considered the slepton tranvcrsc momentum 
distribution and resummed the final state gluons. In- 
cluding these results we have performed a detailed 
comparison with the Herwig and Susygen/Pythia 
Monte Carlo generators. Significant discrepancies re- 
main, similar to the case of SM Z-boson production. 
This needs to be remedied. Finally including the slep- 
ton decay branching ratios we determine the like-sign 
dimuon event rates at the LHC and Tevatron for var- 
ious SUSY spectra. This should enable a detailed in- 
vestigation by the experimental groups. 
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TABLE IV: Number of like-sign dimuon events, fJ, n + X, from cascade decays of fi and at the LHC (Tevatron) with 
an integrated luminosity of 10fb _1 (lfb -1 ) for the six mSUGRA SPS scenarios. The expected number of events include 
also the charge conjugated final states. Events which are kinematically forbidden or with BR < 0.5% or where there are 
less then 0.1 events are denoted by a dash. 
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APPENDIX A: HIGH-p T SLEPTON 
DISTRIBUTION AT 0(a s ) 

The transverse momentum spectrum in the high px 
region, i.e. pt > to, can be calculated perturbatively. 
According to [70, [90] the perturbative hadron level 
cross section for the parton level process 



o(pi) +b(p 2 ) 
can be written as: 

da pert 



L(q) + X 



(A.l) 



dp T dy 

2p T 
S 



1 dx\ 



da 



x* A X\ X-± 
1 



dt 



dx 2 , / »w / \ a do 

+ JaA{X 1 )f b B{X 2 )S— r {Xi = Xi) 

x' B X 2 -xj dt 

(A.2) 

Here, parton a with momentum p\ = X\P\, and par- 
ton b with momentum p 2 = x 2 P 2 produce a slepton 
(sneutrino) (L) with momentum q and mass to and 
an additional particle, X. Parton a originates from 
hadron A with momentum Pi and parton b originates 
from hadron B with momentum P 2 . f a A{x\) denotes 
the parton density function for the parton a in hadron 



A with momentum fraction x\. The other hadron level 
variables as functions of pr are 

S=(Pi+ P 2 f, 

X * A/B = e ±^x/7+, x° 1/2 = e^^rV/S, 

xt /2 = e ±y ^{m?+pl)/S, 

x i/2 = ( x 2/i x t/ 2 " ™?/S)/{x 2/l - x+ fl ) . (A.3) 

Here y is the slepton rapidity and x^ 2 are the parton 
momentum fractions in the limit pr/m — * 0. The 
partonic variables arc given by 





s = 


xix 2 S, 










i = 


TO 2 




Xl 

x\\ 




TO 2 


x* 2 ) 


u = 


TO 2 


(- 


x 2 
x°\ 




„ fi. 

TO 2 



(A.4) 



The transverse momentum spectrum of the slepton 
(sneutrino) is at 0(a s ) due to the gluon radiation and 
Compton-like processes. For gluon radiation, the rel- 
evant parton level processes arc 



q{Pi) + <l'(P2) -> L(q) +g, 
q'(Pi) + q{P2) -> L(q) + g, 



(A-5) 



where q and q' are not necessarily the same flavour. 
The resulting differential partonic cross section is 
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given by 



We obtain the differential partonic cross sections: 



di 



12 



dt 
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h2 — - 
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(A.6) 



and s, i, u are the parton level Mandelstam variables 
(|A.4[) . The parton level Compton-like processes with 
a = q,q' , b = q' ,q and g denoting a gluon are: 



ff(Pi) + a{P2) -> £(?) + &• 



(A.7) 



.da" 3 


a s X' 2 
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- - H 


■»( 


/ fn 2 u 


s — — = 

dt 


48s 
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v. St 




a s A' 2 
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u 

\- - - 


f 2 




S = 

dt 


48s 




s 





(A- 



These can then be inserted in Eq. (|A.2[) together with 
the parton density functions to obtain the correspond- 
ing hadronic differential cross sections. 
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